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ABSTRACT 

The X-ray emission of many active galactic nuclei (AGNs) is absorbed between 
15 and 17 A by means of unresolved (inner-shell) transition arrays (UTAs) of 
Fe M-shell ions. The outflow velocities implied by the Doppler shifts of the 
individual UTAs in the spectrum have never before been measured. Thus, the Fe- 
M absorber has been commonly assumed to be part of the ionized AGN outflow, 
whose velocities are readily obtained from more easily measured spectral lines. 
The best spectrum of Fe-M absorption is available from the integrated 900 ks 
Chandra HETGS observations of NGC 3783, in which some Fe-M ions are clearly 
resolved. We measure the velocities of the individual Fe-M ions in NGC 3783 
for the first time. Surprisingly, we find that the Fe-M absorber, most noticeably 
Fe"*"®, Fe"*"^, and Fe"*"^", is not outfiowing at the same velocity as the previously 
known wind. In fact, it appears to be stationary and therefore not part of the 
outfiow at all. It could, alternatively, be ascribed to the skin of the dusty torus. 
This reduces appreciably the mass loss rate estimated for the NGC 3783 outfiow 
and perhaps for other similar sources as well, in which the various Fe-M ions are 
not resolved. 

Subject headings: galaxies: active — galaxies: individual (NGC 3783) — tech- 
niques: spectroscopic — X-rays: galaxies — line: formation 



1. INTRODUCTION 

The X-ray spectrum of many active galactic nuclei (AGNs) viewed directly towards the 
central source (e.g., Seyfert 1 galaxies) shows the continuum fiux absorbed by numerous 
absorption lines. These lines are generally shifted slightly to shorter wavelengths indicating 
an outfiowing wind. The fiux from the central source is believed to drive the wind and 
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to ionize it to the high charge states observed. The X-ray band is uniquely compact as it 
comprises absorption hnes from the full range of charge states from neutral atoms to H-like 
ions of all elements with Z > 6, which helps to constrain the column density distribution 
over a wide range in ionization. The less ionized species absorb by means of inner-shell 
photoexcitation processes, which play an important role in the formation of discrete features 
in the spectrum. 

In particular, the n = 2 to 3 absorption lines (n being the principal quantum number of 
the active electron) by low charge states of Fe, neutral through Fe+^^, fall in the soft X-ray 
band. Each ion can produce dozens of overlapping lines forming unresolved transition arrays 
(UTAs). Sako et al. (2001) reported the first astrophysical identification of these features 
in the spectrum of the luminous quasar IRAS 13349-1-2438 obtained with the Reflection 
Grating Spectrometer (RGS) on board XMM-Newton. The broad absorption trough observed 
between 16-17 A was attributed to a blend of hnes from Fe+^ to Fe"^^^. The basic atomic 
data needed for spectral modeling of the lines was subsequently given in Bchar ct al. (2001). 
Since its initial discovery, copious detections of Fe M-shell absorption have been claimed in 
other soft X-ray sources (e.g.. Pounds et al. 2004; Kaspi et al. 2004; Gallo et al. 2004, to 
name a few recent ones). Since most instruments are not able to resolve the various ionic 
UTAs, for simplicity, the observed, blended feature was dubbed in the literature as the Fe-M 
UTA. 

Perhaps the most conspicuous astrophysical Fe-M UTA is the one found in NGC 3783. 
Both Chandra and XMM-Newton have observed NGC 3783 for very long spectroscopic ex- 
posures (Kaspi et al. 2002; Behar et al. 2003, respectively). Several detailed global models 
have been constructed to explain the rich absorption spectrum of NGC 3783 (Kaspi et al. 
2001; Blustin et al. 2002; Krongold et al. 2003; Netzer et al. 2003). Owing to its ionization 
sensitivity, the Fe-M UTA could have played a major role in constraining these models. How- 
ever, it was realized that the standard photoionization balance models fail tremendously in 
their prediction of the Fe charge state distribution (Netzer et al. 2003; Netzer 2004; Kraemer 
et al. 2004). Therefore, the current models could not produce the Fe-M UTA spectral fea- 
tures consistently with those of other elements. Recently, Krongold et al. (2005) have used 
heavy rebinning of the Chandra spectrum of NGC 3783 to show that the overall position of 
the absorption trough in the Fe-M UTA region appears to shift between low and high states 
of the ionizing continuum. In none of these works were the individual ionic features studied 
closely. 

In this paper, we revisit the integrated 900 ks observation of NGC 3783 obtained with 
the HETGS spectrometer on board Chandra over several months. The same data were used 
before by Kaspi et al. (2002), Netzer et al. (2003), and by Krongold et al. (2003). Presently, 



-3- 



HETGS is the only instrument that can actually resolve the individual ionic features inside 
the Fe-M UTA. Conversely, the spectral resolution of the RGS on board XMM-Newton is 
just insufficient (Sako et al. 2001; Behar et al. 2003). Curiously, the details of the UTA in 
the HETGS spectrum have yet to be studied in detail, perhaps because of the difficulties 
mentioned above in using the Fe-M UTA in global models. The sole goal of this work is 
to use the HETGS spectrum to identify the individual ionic features and to measure their 
velocity with respect to the systemic velocity of NGC 3783. The results obtained are quite 
surprising. 



2. HETGS SPECTRUM OF NGC 3783 

The spectrum used in the present work is the exact same one reduced in Kaspi et al. 
(2002) and shown in their Fig. 1. It consists of five 170 ks HETGS observations of NGC 3783 
obtained during the period 2001 February- June as well as an early 56 ks observation from 
2000 January, all of which are co-added. The -|-lst and —1st diffraction orders of both the 
MEG and HEG grating assembhes are included. The spectra are fluxed and corrected for 
Galactic absorption. The wavelength scale is de-redshifted to the rest frame of NGC 3783 
using z = 0.00976. For more details on the data reduction see Kaspi ct al. (2002). For the 
obvious reason of improving the statistics, we use here the total 900 ks data set. Notwith- 
standing, wc have verified that all of the prominent absorption features, which wc discuss 
below, appear in each short exposure. Moreover, to the extent that it is possible to determine 
with the low statistics of each short exposure, we have found that the positions of the narrow 
absorption features, which are of interest for this work, do not vary between observations. 



3. SPECTRAL MODEL 

The heavy line blending in the soft X-ray band where the outfiow absorbs most sig- 
nificantly makes the direct measurement of individual line positions very difficult. For a 
consistent measurement, therefore, one needs to use an ad-hoc spectral model, which in- 
cludes simultaneously many charge states and all of the absorption lines pertaining to these 
ions in the soft X-ray band. Interested solely in the absorption features of the Fe-M UTA, 
we have constructed a simple model to explain the observed spectrum in the limited range 
of 14.9-17.5 A. The model includes a power-law continuum, all of the Fe-M UTA ions, Fe"*"^^, 
and the two 0-K ions: O"^^ and O"^^, which also have lines in this range. For the Fc ions, 
the model includes all of the n = 2 to 3 spectral lines for absorption from the ground level 
of each ion. The present Fe-M model includes the full set of absorption lines from Behar 
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et al. (2001) and not only the abbreviated set given in Table 1 of that paper. For the 0-K 
ions, the model includes all of the n = 1 to n' resonance transitions with n' < 15000 as well 
as the photoelectric edges, from the ground level. The present model is similar in spirit to 
those used in Sako et al. (2001) and in Behar et al. (2003). 

For the continuum, a power law was fit to the relatively featureless region of the spectrum 
between 2-6 A. The best-fit parameters obtained arc a photon index (F) of 1.45 and a 
normalization at 1 keV of 0.00936 photons s~^ cm"^ keV^^, which are in good agreement 
with previous fits for the low-state spectrum that dominates the present data set (e.g.. Fig. 8 
in Netzer et al. 2003). 

The optical depth Tij{v) due to an absorption line {i — > j) is obtained by: 

Tij{v) ^ Niondijiy) (1) 

where (Tij{i') is the calculated cross section for photoexcitation (in cm^) from i to j and Nion is 
the ionic column density towards the source (in cm~^) for which we fit. The photoexcitation 
cross section is given by: 

<^iA^) = zrzfiM^) (2) 

Here, e and rrie are the electron charge and mass, c is the speed of light, fij is the line os- 
cillator strength, which wc have calculated with the Hebrew University Lawrence Livcrmorc 
Atomic Code (HULLAC, Bar-Shalom et al. 2001). Finally, (j){iy) is the normalized line profile 
described here by a Voigt function that is the convolution of the (Gaussian) Doppler broad- 
ening due to temperature and to turbulence, which we fit for, and the (Lorentzian) natural 
width of the line, which we have calculated with HULLAC. In observed AGN outflows, tur- 
bulence broadening {vturb > 100 km s~^) totally dominates over temperature broadening. 
Here, we use Vturb = 170 km s~^. Generally, in the natural widths of inner-shell excited lines, 
the autoionization decay rate {A"- ^ 10^'* s~^) dominates over the radiative rate. Neverthe- 
less, our calculations include both. HULLAC is a multi-configuration, relativistic code most 
suitable for the inner-shell excitations of the complex atoms of interest in this study. For 
more details on the Fe-M atomic calculations see Behar et al. (2001). 
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4. DETERMINATION OF OUTFLOW VELOCITIES 

The mean velocity for the NGC 3783 outflow was found by Kaspi et ah (2002) to be 
—590 ± 150 km s^^. Consequently, our flrst attempt to flt the absorption in the HETGS 
spectrum of NGC 3783 in the vicinity of the Fe-M UTA was carried out using a uniform 
outflow velocity of -590 km s"^. The results are shown in Fig. 1. Indeed, the flt for 0+^ 
(He 5 and He e most prominent at 17.40 and 17.20 A), 0+^ (Ly (3 and Ly 7 at 16.01 and 
15.18 A), and Fe"*"^^ (at 15.01 and 15.26 A) is seen to be very good. On the other hand, 
the velocity of -590 km s~^ is clearly inappropriate for the prominent Fe-M troughs. The 
Fe-M features in the model are seen to be systematically shifted towards shorter wavelengths 
compared with their observed positions, most evidently at 16.50, 16.33, and 16.15 A. 

The above result is already a strong indication that the Fe-M absorber is moving at a 
much slower velocity than the more highly ionized gas, if it is moving at all. In Table 1 we 
list the identiflable discrete absorption features in the spectrum. For the mere objective of 
obtaining centroids in order to get a rough idea of the velocities, we have fltted each feature 
with a local continuum and a Gaussian-shaped trough employing the same method used for 
the wavelength determination in Kaspi et al. (2002). It can be seen that the positions of the 
best observed Fe-M ions are consistent with outflow velocities that are signiflcantly lower 
than — 590 km s~^. The velocities in the table have to be treated with caution for two main 
reasons: (i) The absorption line proflle in many cases is not a simple Gaussian, but rather a 
complex distribution of narrow velocity components (see also Kaspi et al. 2002; Gabel et al. 
2003) (ii) The Fe-M wavelengths have not been conflrmed in the laboratory as discussed in 
detail in §5. The spectral resolution of the HETGS is high enough so that blending among 
Fe-M ions is negligible as can be seen by comparing the 3rd and 4th columns in Table 1, 
which show that the individual-ion Fe-M wavelengths and the full-model wavelengths agree 
to within a few mA. 
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Table 1. Best-fit velocities and column densities for ions detected in the 14.9-17.5 A 
region of the HETGS spectrum of NGC 3783. 



Ion 


'^Observed 




'^Model 


Outflow Velocity 


Ion Column Density 




(A) 


(A) 


(A) 


(km s~^) 


(1016 cm-2) 




15.144 ± 0.003 <^ 


15.176 


15.146 


-632 ± 59 


400 ± 60 




15.970 ± 0.005 ^ 


16.006 


15.973 


-665 ± 94 




0+6 


17.351 ± 0.005 


17.395 


17.361 


-770 ± 81 


110 ± 20 




17.161 ± 0.005 


17.199 


17.165 


-675 ± 89 






17.048 ± 0.005 


17.084 


17.050 


-630± 88 




Fe+16 


14.980 ± 0.003 


15.013 


14.985 


-659 ± 60 


3.0 ± 0.5 




15.231 ± 0.002 ^ 


15.261 


15.234 


-590 ± 39 




Fe+15 


15.231 ± 0.002 ^ 


15.250 


15.234 


-374 ± 39 


0.6 ± 0.2 


Fe+14 


15.322 ± 0.006 


15.316 


15.317 


+ 118 ± 118 


0.3 ± 0.1 


re 


lo.ooy zt u.uuo " 


iO.OoU 


io.ooU 




^ A 4- n A 


Fe+12 


15.844 ± 0.014 s 


15.848 


15.848 


-76 ± 265 


1.0 ± 0.3 


Fe+ii 


15.970 ± 0.005 ^ 


15.967 


15.973 


+56 ± 94 


3.0 ± 2.0 


Fe+io 


16.154 ± 0.005 


16.150 


16.150 


+74 ± 93 


4.0 ± 0.7 


Fe+9 


16.329 ± 0.005 


16.339 


16.340 


-184 ± 92 


5.5 ± 0.7 


Fe+8 


16.496 ± 0.004 


16.510 


16.509 


-252 ± 73 


4.0 ± 0.5 


Fe+7 h 




16.655 






3.0 ± 1.0 


Fe+6 h 




17.095 






2.0 ± 0.7 


Fe+5 h 




17.218 






1.5 ± 0.6 


Fe+4 h 




17.291 






1.0 ± 0.6 


Fe+3 h 




17.405 






< 0.8 



"^For Fe-M ions, these are centroids of the deepest feature in each ionic spectrum. 
''Centroid in full two- velocity model (Fig. 2). 

'^Estimated by (Ao6ser-yed ~ ^Rest)c/^Rest- Errors reflect 90% confidence intervals. 
•^Blend of 0+'^ and Fe+^^ at 15.19 A fitted for two distinct features. 

''Unresolved blend of 0+^ and Fe"*"^^ reflected in the large error on the Fe"*"^^ column density. 
''Unresolved blend of Fe+^^ and Fe+^^ 



- 7- 



^arge uncertainty in centroid identification. See spectrum in Figs. 1 and 2. 

^Indirect identification based on best-fit model; Strongly blended with the high-n lines of 0+^. 
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In order to further illustrate the low velocities associated with the Fe-M absorber, we 
have fitted the HETGS spectrum of NGC 3783, but now the velocity of the Fe-M ions was 
allowed to be different from that of the other charge states. For simplicity, however, the 
velocities of all of the Fe-M ions were tied together. The best fit was obtained by setting the 
Fe-M velocity to zero as shown in Fig. 2. The individual contributions of Fe^'' to Fe+^^ to 
the total model are plotted in Fig. 3. Despite some blending, individual ions can be easily 
discerned in the spectrum. As can be seen, centroids of the strongest absorption features in 
the model are determined almost exclusively by a single UTA, which reinforces the validity 
of the direct velocity estimates in Table 1. The most conspicuous Fe-M ions in the data are 
Fe"*"^, Fe"*"^, and Fe"'"^°. These ions in particular produce UTAs that are only very weakly 
blended with O lines or with other Fc ions as demonstrated in Fig. 3. Particularly for these 
ions, the fit of the model with zero outfiow velocity can be seen in Fig. 2 to be very good. The 
obvious velocity discrepancies for these ions in Fig. 1 have now disappeared. Conclusions 
regarding higher and lower Fe-M charge states arc not as obvious, although the data still 
favor the stationary scenario over the — 590 km s^^ outfiow for these ions as well, as was 
already indicated by the direct velocity estimates presented in Table 1. 



5. DISCUSSION 

The outfiow of NGC 3783 is known from high-resolution UV observations to have four 
resolved velocity components (Gabel et al. 2003). In the X-rays, the four components are 
unresolved. However, the mean X-ray absorption line profile is generally in agreement with 
the UV profile, as demonstrated in Fig. 2 of Gabel et al. (2003). The lowest velocity UV 
absorber reported in Gabel et al. (2003) is at —548 km s~^, where the absorption troughs 
rise sharply. On the other hand, from our own inspection of the HST and FUSE spectra 
of NGC 3783 we find a much slower absorber, the most prominent lines of which are Lya 
and Ljl3 observed at an infall velocity of about +100 km s~^ (see also Fig. 1 in Gabel et al. 
2003). As far as we know, this together with the present X-ray results represents the first 
discovery of a high column density, almost stationary X-ray component in any AGN. In 
the intensive XMM-Newton study of NGC 3783, the low charge states of oxygen (neutral O 
and O"*"^) have actually been detected at zero velocity as well (Behar et al. 2003, Table 1 
therein). However, both lines were stated to be blends and in any event neither feature was 
statistically significant. The present detection is much more unambiguous. 

The stationary or low- velocity component could possibly be an extended, slightly ionized 
part of the interstellar medium (ISM) of NGC 3783. It might be absorption created by the 
extended narrow fine region of NGC 3783. It could also be gas evaporating off of the molecular 
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torus. A scenario in which radiation pressure is inefficient has been suggested for NGC 3783 
by Chelouche & Netzer (2005). The outflow velocity has obvious consequences on the mass 
outflow rate, which for NGC 3783 has been estimated to exceed a solar mass per year (Behar 
et al. 2003). In fact, Netzer ct al. (2003) have estimated the mass loss rate in this target due 
to the low ionization component alone to be as high as 75 solar masses per year and scaling 
linearly with the radial filling factor of the fiow, which could be as low as 0.1 (Netzer et al. 
2003). All of these estimates for the mass loss rate are based on the high velocity of the 
outflow (~ —590 km s~^). If indeed, much of this gas is not moving as fast, it would reduce 
the estimated mass loss rate considerably. 

Outflow velocities for individual Fe-M ions have not been carried out for other AGN tar- 
gets, mostly because the available observations generally do not resolve the individual ions in 
the spectrum, even when the blended Fe-M UTA feature is very clear (e.g., Sako et al. 2001). 
The present result raises the suspicion that perhaps in other targets as well, this component 
of the absorber is not part of the AGN outflow. Of course, this will need to be checked in 
the future for each target separately. However, if indeed this is the case, the ionization state 
diagnosed for these absorbers based on the position of the UTA is underestimated, as the 
high velocities assumed for these targets are being confused with higher charge states that 
are actually moving at lower velocities. 

Before concluding, we need to mention one disturbing caveat about the present measure- 
ment. The rest frame wavelengths of the resonant n = 2 to 3 absorption lines, which we use 
to determine the Fc-M outflow velocity, arc taken from theoretical calculations (Behar ct al. 
2001). Measuring these wavelengths in the laboratory is rather challenging and indeed, they 
have not yet been confirmed. As far as we know, Fe-M UTA absorption has been observed 
in the laboratory once by Chenais-Popovics et al. (2000). However, the spectral resolution 
of that measurement was not sufficiently high to pin down the positions of individual charge 
states. 

The errors on the present velocity estimates stem directly from the position uncertainties 
of the strongest fines in each UTA. At 16 A, an uncertainty in the rest frame wavelength 
of 30 mA translates into an uncertainty of 560 km s~^ in the outfiow velocity. Indeed, the 
scatter among wavelength calculations by different state-of-the-art atomic codes, which is 
one way to get a rough idea of their expected uncertainty, can reach and even slightly exceed 
30 mA. This, e.g., is the case for the inner-shell Ka absorption lines of oxygen (Behar & 
Kahn 2002; Schmidt et al. 2004). However, even if there are such uncertainties associated 
with the Fe-M rest frame wavelengths that we use, there is no particular reason for them 
to all be off in the same direction as to indicate a uniformly erroneous velocity. Quite the 
contrary, with the parametric potential method used by HULLAC (Bar-Shalom et al. 2001), 
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the atomic potential is calculated ab initio and separately for each ion, so that errors can be 
expected to be unconnected. Indeed, a comparison between EBIT laboratory measurements 
and HULLAC wavelengths for the strongest line (the same 2p - 3d transition that dominates 
the Fe-M UTAs) in each of the Fe-L ions reveals that Ahullac ~ Aebit = 0, —5, +6, —12, 
+8, +12, +27, and +2 mA, respectively, for Fe+^^ - Fe+^s (Brown et al. 1998, 2002). Note 
that the laboratory wavelengths themselves are associated with errors of a few mA. In short, 
the HULLAC errors are generally much smaller than 30 mA and are certainly not uniform. 
Having said that, until the Fe-M UTA wavelengths are measured directly, this possibihty can 
not be rigorously ruled out. If it turns out that after all, the apparent discrepancy between 
the outflow velocities of the Fe-M low charge states and faster components of the outflow 
is due solely to the inaccurate Fc-M atomic data we have used, then the HETGS data and 
Table 1 could serve as (the best currently available) wavelength calibration for the Fe-M 
wavelengths. 

6. CONCLUSIONS 

We have measured the velocity of individual Fe-M ions in the X-ray outflow of NGC 3783 
for the first time. This measurement has never been carried out before for NGC 3783 nor for 
any other target. Based on the theoretical rest frame wavelengths available to us, we find 
that the Fe-M ions arc not moving as fast as the more ionized Fe species, such as Fe^-'^^, or 
as other highly charged species, which comprise the high-velocity outfiow of NGC 3783. In 
fact, the Fe-M absorber may not be moving at all along our line of sight. We speculate that 
this stationary gas may pertain to the ISM of NGC 3783, to its narrow line region, or it may 
be gas evaporating off of the skin of the molecular torus. In order to examine whether this is 
a distinctive quality of NGC 3783 or if there is a fundamental problem with the atomic data, 
more high-quality spectra will need to be obtained for similar targets, which show deep, 
resolved Fe-M UTA troughs. 
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Fig. 1.— HETGS spectrum of NGC 3783 in the Fe-M UTA region fitted with a uniform- 
velocity (-590 km s~^) outflow model. The model fits some of the data very well, but not 
the Fe-M UTA features, most noticeably Fe"*"*, Fe"*"^, and Fe"*"^*^ for which this velocity is 
clearly too high. 
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Fig. 2.— HETGS spectrum of NGC 3783 in the Fe-M UTA region fitted witli a model in 
wliicfi O"^^, O^^, and Fe"*"^^ are outflowing at a velocity of —590 km s~^, while the Fe-M ions 
are at rest with respect to the host galaxy. This model fits the observed Fe-M UTA features 
appreciably better than the uniform- velocity model shown in Fig. 1. 
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Fig. 3. — Individual ionic contributions to the full best-fit model plotted in Fig. 2. For each 
ion, the strongest absorption feature can be easily discerned and its centroid position in the 
full model (red curve) is essentially unaffected by blends. 



